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We present the first observation of the decay B — » 
J/ip(f>K. Using 9.6 x 10 6 BB meson pairs collected with 
the CLEO detector, we have observed 10 fully reconstructed 
B —* J /tp (j> K candidates, whereas the estimated background 
is 0.5 ± 0.2 events. We obtain a branching fraction of B(B — > 
J/tp<t>K) = (8.8t^[stat] ± 1.3[syst]) x 10~ 5 . This is the first 
observed B meson decay requiring the creation of an addi- 
tional ss quark pair. 



An observation of a B meson decay requiring the cre- 
ation of an additional ss quark pair in the final state 
would enhance our understanding of strong interactions 
in the final states of B decays. Previous studies of 
such processes involved searches for the "lower vertex" 
B — > DfX transitions however no signal was ob- 
served. The decay B — > J/ijj <p K Q can occur only if an 
additional ss quark pair is created in the decay chain be- 
sides the quarks produced in the weak b — > ccs transition. 
The B — > J/tjxpK transition most likely proceeds as a 
three-body decay (Fig. [j]) . Another possibility is that the 
B — > J/tp (f> K decay proceeds as a quasi-two-body decay 
in which the J / tp and <f> mesons are daughters of a hybrid 
charmonium state ||. 

We searched for B+ -> J/tp 4> K+ and B° -> J/tp 4> K° s 
decays, reconstructing J/tp — » £ + £~, <p — ► K + K~, and 
Kg — > 7r + 7r _ . Both e + e~ and modes were used for 

the J/tp reconstruction. The data were collected at the 
Cornell Electron Storage Ring (CESR) with two config- 
urations of the CLEO detector, called CLEO II |J] and 
CLEO II. V. The components of the CLEO detector most 
relevant to this analysis are the charged particle tracking 
system, the Csl electromagnetic calorimeter, the time-of- 
flight system, and the muon chambers. In CLEO II, the 
momenta of charged particles are measured in a track- 
ing system consisting of a 6-layer straw tube chamber, 
10-layer precision drift chamber, and 51-layer main drift 
chamber, all operating inside a 1.5 T solcnoidal magnet. 
The main drift chamber also provides a measurement of 
the specific ionization loss, dE/dx, used for particle iden- 
tification. For CLEO II. V, the innermost wire chamber 
was replaced with a three-layer silicon vertex detector || . 
Muon identification system consists of proportional coun- 
ters placed at various depths in the steel absorber. 
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The results of this search are based upon an integrated 
luminosity of 9.1 fb _1 of e + e _ data taken at the T(45) 
energy and 4.4 fb _1 recorded 60 MeV below the T(45) 
energy. The simulated event samples used in this anal- 
ysis were generated using GEANT-based Q simulation 
of the CLEO detector response. Simulated events were 
processed in a similar manner as the data. 

When making requirements on such kinematic vari- 
ables as invariant mass or energy, we took advantage of 
well-understood track and photon-shower covariance ma- 
trices to calculate the expected resolution for each com- 
bination. Therefore we extensively used normalized vari- 
ables, which allowed uniform candidate selection criteria 
to be used for the data collected with the CLEO II and 
CLEO II. V detector configurations. 

The normalized invariant mass distributions for the 
J/tp — » i + £~ signal in data are shown in Fig. ||. We re- 
quired the normalized invariant mass to be from —10 to 
+3 (from -4 to +3) for the J/tp e + e-{J/tp -» ^ + ^) 
candidates. The resolution in the £ + £~ invariant mass 
is about 10 MeV/c 2 . To improve the energy and mo- 
mentum resolution of a J/tp candidate, we performed a 
fit constraining the mass of each J/tp candidate to the 
world average value [Q. 

Electron candidates were identified based on the ra- 
tio of the track momentum to the associated shower en- 
ergy in the Csl calorimeter and specific ionization loss 
in the drift chamber. The internal bremsstrahlung in 
the J/tp — ► e + e~ decay as well as the bremsstrahlung 
in the detector material produce a long radiative tail in 
the e + e~ invariant mass distribution and impede effi- 
cient J/tp — > e + e~ detection. We recovered some of the 
bremsstrahlung photons by selecting the photon shower 
with the smallest opening angle with respect to the di- 
rection of the e ± track evaluated at the interaction point, 
and then requiring this opening angle to be smaller than 
5°. The addition of the bremsstrahlung photons re- 
sulted in a relative increase of approximately 25% in the 
J/tp — > e + e~ reconstruction efficiency without adding 
more background. 

For the J/tp — > reconstruction, one of the muon 

candidates was required to penetrate the steel absorber 
to a depth greater than 3 nuclear interaction lengths. 
We relaxed the absorber penetration requirement for the 
second muon candidate if it was not expected to reach a 
muon chamber either because its energy was too low or 
because it pointed to a region of the detector not covered 
by the muon chambers. For these muon candidates we 
required the ionization signature in the Csl calorimeter 
to be consistent with that of a muon. Muons typically 
leave a narrow trail of ionization and deposit approx- 
imately 200 MeV of energy in the crystal calorimeter. 
Hadrons, on the other hand, quite often undergo a nu- 
clear interaction in the Csl crystals that have a depth of 
80% of a nuclear interaction length. Compared to impos- 
ing the absorber penetration requirement on both muon 



2 



candidates, this procedure increased the J/ip — > 
reconstruction efficiency by 20% with 80% increase of 
background. 

We required that the charged kaon candidates have 
dE/dx and, if available, time-of-flight measurements that 
lie within 3 standard deviations of the expected values. 

If for the B — > J/ip <p K decays we assume a uniform 
Dalitz distribution and isotropic decays of J/ip and <p 
mesons, then the expected efficiency of the combined 
dE/dx and time-of-flight selection is approximately 90% 
per kaon candidate. The dE/dx measurements alone pro- 
vide the K/ir separation of more than 4 standard devia- 
tions for 92% of the (p daughter kaons and for 64% of the 
"bachelor" kaons from B decay. We selected (p — » K + K~ 
candidates by requiring the K + K~ invariant mass to 
be within 10 MeV/c 2 of the cp mass [@. We did not 
use the normalized K + K~ invariant mass because the 
mass resolution (1.2 MeV/c 2 ) is smaller than the (p width 
(4.4 MeV) §. 

The Kg candidates were selected from pairs of tracks 
forming well-measured displaced vertices. The resolution 
in 7r + 7r~ invariant mass is about 4 MeV/c 2 . We required 
the absolute value of the normalized tt + tt~ invariant mass 
to be less than 4, then we performed a fit constraining the 
mass of each Kg candidate to the world average value . 

The B — > J/ip (p K candidates were selected by means 
of two observables. The first observable is the difference 
between the energy of the B candidate and the beam 
energy AE = E(J/ip) + E{4>) + E{K) - E heam . The 
resolution in AE for the B — > J/ip (p K candidates is ap- 
proximately 6 MeV. The second observable is the beam- 
constrained B mass M(B) = */ E^ — p 2 (B), where 
p(B) is the absolute value of the B candidate momentum. 
The resolution in M(B) for the B — > J /ip <pK candidates 
is about 2.7 MeV/c 2 ; it is dominated by the beam energy 
spread. The distributions of the AE vs M(B) for B + — > 
J/ip(pK+ and B° -> J/ipcpK s are shown in Fig. |. We 
used the normalized AE and M(B) variables to select the 
B — » J/ip(pK candidates and defined the signal region 
as \AE/a(AE)\ < 3 and \(M(B) - M B )/a(M{B))\ < 3. 
We observed 8(2) events in the signal region for the 
B+ -> J/iP(PK+ (B° -> J/ipcj)K s ) mode. Consider- 
ing that K° can decay as K s or as K^, and also taking 
into account B(K S — > 7r + 7r~) and the difference in recon- 
struction efficiencies, we expect to observe on average 4.3 
B + — > J/ip cp K + candidates for every B° — ► J/ip(p K s 
candidate. 

The Dalitz plot and the cosine of helicity angle dis- 
tributions for the 10 B — > J/ip(pK signal candidates 
are shown in Figs. [| and ||. The helicity angle for 
J/ip £ + £~ decay is defined as the angle between a 
lcpton momentum in the J/ip rest frame and the J/ip 
momentum in the B rest frame. An analogous defini- 
tion was used for the <p — > K + K~ decay. No conclu- 
sion can be drawn yet either about the J/ip and the <p 
polarizations or about the resonant substructure of the 



B — > J /ili (pK decay. If the J / ip and (p mesons are the 
products of the hybrid charmonium ip g decay, then the 
J/ip<p invariant mass is expected to be below the DD** 
threshold (4.3 GeV/c 2 ) because ip g — > DD** decay is 
likely to dominate above the threshold || . The J/ip (p in- 
variant mass is above 4.3 GeV/c 2 for all 10 B — > J/ip (pK 
candidates thus difavoring the hybrid charmonium dom- 
inance scenario. 

The background can be divided into two categories. 
First category is the combinatorial background from 
T(45) — > BB and continuum non-BB events. Sec- 
ond category is the background from non-resonant B — > 
J/ip K+ K- K decays. 

The combinatorial background from T(45) — * BB 
events was estimated using a sample of simulated events 
approximately 32 times the data sample; events contain- 
ing aJ3-» J/ip K + K~ K decay were excluded. We esti- 
mated the background from T(45) — > BB decays to be 
0.25^Q;pg events. In addition, we specifically considered 
B -> J/ip K* 7T+ with K* -> Kir~ and B -> J/ip p° K 
decays because the beam-constrained B mass distribu- 
tion for these modes is the same as for the B — > J/ip (pK 
decays. Using data and simulated events, we verified 
that those backgrounds are rendered negligible by the 
kaon identification, <p mass, and AE requirements. The 
combinatorial background from the continuum non-BB 
events was estimated using simulated events and the data 
collected below BB threshold. We found the continuum 
background to be negligible. 

To estimate the background contribution from the non- 
resonant B — * J/ip K + K~ K decays, we reconstructed 
B+ -> J/iJj K+ K- K+ and B° ' -> J/iP K+ R- K° s 
candidates in data requiring \M(K + K~) — M$\ > 
20 MeV/c 2 to exclude B — > J/ip (pK events. We observed 
7 B — > J/ip K + K~ K candidates with the estimated BB 
combinatorial background of 2.8 events. We estimated 
the mean background from B — > J/ip K + K~ K decays 
for the B — * J/ip<pK signal to be 0. 27to'i7 events; we 
assumed that B — > J/ip K + K~ K decays according to 
phase space. 

In summary, the estimated total background for the 
combined B — > J/ipcpK signal is 0.52^p' 2 g events. 

We evaluated the reconstruction efficiency using a sam- 
ple of simulated B — > J/ip (p K decays. We assumed a 
uniform Dalitz distribution and isotropic decays of J/ip 
and <p mesons; these assumptions are consistent with data 
(Figs. |J and ||). The reconstruction efficiency, which does 
not include branching fractions of daughter particle de- 
cays, is (15.5 ±0.2)% for the B+ -> J/ipcpK+ mode and 
(10.3 ± 0.2)% for the B° -> J/ip(pK° s mode. The recon- 
struction efficiency is close to zero at the edges of phase 
space where either (p or K meson is produced nearly at 
rest in the laboratory frame. Thus, the overall detec- 
tion efficiency would be much smaller than the above 
values if the B — > J/ip (pK decay is dominated by ei- 
ther a J/ip K resonance with a mass around 4.3 GeV/c 2 
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or a J/ijj<f> resonance with a mass around 4.8 GeV/c 2 . 
No such resonances are expected. To assign the system- 
atic uncertainty due to the decay model dependence of 
the reconstruction efficiency, we generated two additional 
samples of simulated B — > J/ip(j>K events. One sam- 
ple was generated with a uniform Dalitz distribution for 
B — > J/iptf>K and 100% transverse polarization for J/tp 
and 4>. The other sample was generated assuming the <j) 
and K mesons to be daughters of a hypothetical spin-0 
resonance with mass 1.7 GeV/c 2 and width 100 MeV. We 
estimated the relative systematic uncertainty due to the 
decay model dependence of the reconstruction efficiency 
extraction to be 7%. 

For the branching fraction calculation we assumed 
equal production of B + B~ and B°B pairs at the 
T(4S) resonance and B{B+ -> J/ifjcj>K + ) = B{B° -> 
J/iP(j)K ) = B(B -► J/ip<f>K). We did not assign any 
systematic uncertainty due to these two assumptions. 
We used the world average values of B{Jhb — > 
B(<f> -» K+K~), and B(K° S -» tt+tt") @. We used 
the tables in Ref. M to assign the 68.27% C.L. inter- 
vals for the Poisson signal mean given the total num- 
ber of events observed and the known mean background. 
The resulting branching fraction is B(B — > J/tptpK) — 
(8.8t^[stat] ± 1.3[syst]) x 10~ 5 . 

The systematic error includes the uncertainty in the 
reconstruction efficiency due to decay modeling plus the 
uncertainties in track finding, track fitting, lepton and 
charged-kaon identification, Kg finding, background sub- 
traction, uncertainty in the number of BB pairs used for 
this measurement, statistics of the simulated event sam- 
ples, and the uncertainties on the daughter branching 
fractions B(J/ip -> and B(<j> -» K+K~) §. We 

estimated the total relative systematic uncertainty of the 
B(B — > J/ip(f)K) measurement to be 15%. 

In conclusion, we have fully reconstructed 10 B — > 
J /ip (j>K candidates with a total estimated background 
of 0.5 events. Assuming equal production of B + B~ 

and B°B pairs at the T(45) resonance and B(B + —> 
J/ip<pK+) = B(B° -> J/^0A°) = B(B -> J/tp<pK), 
we have measured B(B -> J/ip<f>K) = (8.8^0 [stat] ± 
1.3(syst)) x 10~ 5 . This is the first observed B meson 
decay requiring the creation of an additional ss quark 
pair. 
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FIG. 1. Most likely B — ► J/tp<j)K decay mechanism. 
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FIG. 2. Normalized invariant mass of the (a) J/tp — > e + e~ and (b) J/ip — » candidates in data. The solid line represents 

the data taken at the T(4S) energy; the dashed line represents the luminosity-scaled off resonance data showing the level of 
background from non-SB events. 
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FIG. 3. The AE vs M (B) distribution for (a) B + -> J/tp(j)K + and (c) B° -> J/tp(f)K s candidates in data. The signal 
candidates, selected using normalized AE and M(B) variables, are shown by filled circles. The M(B) distribution for (b) 
B + — » J/ip(f>K + and (d) B° — > J/ip<j)K% candidates satisfying |A_E/ct(A_E)| < 3; the shaded parts of the histograms represent 
signal candidates. 
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FIG. 4. Dalitz plot for the 10 B — » J/ip <j> K candidates in data. The kinematic boundary is represented by solid line. 
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FIG. 5. The distribution of the cosine of helicity angle for Jjip 
the 10 B — > J/rp <j> K candidates in data. 
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